Sheet 5 


Determine if each of the following statements is true or not. If it is not, modify it to become 
true. 

a) When dealing with intensity values as random variable. Its variance is proportional to the 
image contrast. 

b) An intensity mapping on a 8 bit image has the form s = T(z ) = 2 8 - z, where z is the 

pixel intensity in the input image and s is the pixel intensity in the output image, washes 
the image. 

c) In image processing terminology, image registration is the process of aligning input images 

to a master image so that comparison and measurements could be done across all the 
images. 

d) In digital image processing, when dealing with image transform, it is .always possible to 
formulate the transform using matrix multiplications 

e) In digital image processing spatial operations are usually more computationally intensive 
compared to transform operation 


2 . 


3. 


4. 


a) 


Prove that the Fourier transform kernels are separable and symmetric 

r(x,y,u, v) = e~i 2n(r M + ~ffa 

s{x.y,u,v) = - L e K*$h%y 
MN 

Show that the two dimension transform with separable, symmetric kernels can be computed bv 
(1) computing 1-D transform along the individual rows (columns) of the input, followed by (2 
computing 1-D transform along the columns(rows) of the result of step (1) 

An intensity mapping on a 8 bit image has the form s = T(z ) = 2 8 - z where z is the 
eS or 5 '^ the inPUt an - 5 15 the PiXe ' intenSitV in the outpu * irna 8 e ' has the 

I- Washing the image 

II. ^creasing the contrast of the image 

III. Decreasing the contrast of the image 

IV. Getting the negative of the image 

When dealing with intensity values as random variable, positive third moment means 
Pixel values have bias to values higher than the mean 

II- Pixel values have bias to values smaller than the mean 

Whpn hIT ValU ! S h3Ve bi3S t0 be uniform, V distributed around the mean 

I PixEue 'h tenS h V Va ' UeS 35 rand ° m Variab ‘ e ' neg3tiVe th ' rd moment means 
xel values have bias to values higher than the mean 

Pixel values have bias to values smaller than the mean 

Pixel values have bias to be uniformly distributed around the mean 


b) 


d) When dealing with intensity values as random variable, zero third moment means 
Pixel values have bias to values higher than the mean 
"• Plxel va| ues have bias to values smaller than the mean 
III. Pixel values have bias to be uniformly distributed around the mean 
n igita image processing, neighborhood and intensity mapping are 
I- Spatial operations 
II. Transform operations 

HI. Spatial and transform operations 


annl H , T ati ° n * ^ mappm g of the coordinates of each pixel in an input image tc 

differ!' SP c /r0tated '-' ) pixel in the ° utput ima § e - This ™pp,ng can be done in twc 

d ffere! he! T rd ' maPP ' n8 9nd inverse ' ma PP' n g- Explain the two ways in contrasting the 
difference between them. Which way you prefer? Why? 

What is image registration? Why and when they are necessary? 
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Sheet 5 solution 


Item 

answer 

The true statement if False 

A 

T 


B 

F 

An intensity mapping on a 8 bit image has the form s = T(z) = 2 8 — z, 
where z is the pixel intensity in the input image and s is the pixel intensity in 
the output image, performs the negative of the image. 

C 

T 


D 

F 

It is only possible to formulate the transform in terms of matrix multiplication 
if the following two conditions apply 

• Both the forward and the reverse kernels are separable and 
symmetric 

• The operated image is square M by M 

E 

F 

In digital image processing transform operations are usually more 
computationally intensive compared to spatial operation since in most 
transform operation each output pixel value is calculated from all pixels in 
the input image. Also in most cases the process requires to do inverse 
transformations. 


_ o c ux v y ~\ 
r(x, y,u,v) = e J ^ m + n> 

= e - i2n ^e~ i2n ^ 

= r 1 {x,u)r 2 {y,v) 

Hence the kernel is separable 

It's also symmetric (the same form for the two parts) 

^ % ! j?Tr( ux 

The steps can be applied to the reverse kernel: s(x,y,u,v) = — e 1 n } 


3. Using the forward transform and applying the separable property 

M— 1 N-l 


T (u, v) = ^ / (x,y)r(x, y, u, v) ■ 


: X=Qy=0 

M- 1 N-l 


x=0 


= X n X f( x ’ y ) r i( y - 

y=0 

' W-l ; 

= Zk*,fr 2M 


Where 


* 


y=0 
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N - 1 


T(x,v) - ^ f(x, y)r 2 (y,v) 
y=o 

For a fixed value of x, this equation is recognized as the 1-D transform along one row of f (x,y ). By 
letting x vary fromO to M -1 we compute the entire array T(x,v ). Then, by substituting this array 
into the last line of the previous equation we have the 1-D transform along the columns of T(x,v ). 
In other words, when a kernel is separable, we can compute the 1-D transform along the rows of 
the image. Then we compute the 1-D transform along the columns of this intermediate result to 
obtain the final 2-D transform, T (u,v ). We obtain the same result by computing the 1-D transform 
along the columns of f (x,y ) followed by the 1-D transform along the rows of the intermediate 
result. 

This result plays an important role in Chapter 4 when we discuss the 2-D Fourier transform. From 
Eq. (2.6-33), the 2-D Fourier transform is given by 



M— 1 N-l 


t(u, v) = f(. x >y)e~ J2n ^ + ^ 


x=0 y - 0 

Since the transform kernel is separable, we can write 


M-l N-l 




t(u,v) = 2^ 2 l, f( x >y} e ' 

x=0 y = 0 

M—l N-l 

T(u,v ) = ^ ^ /(x,y)e _ 727 r ^AT^ 

x=0 y = 0 

M-l 

T(u,v ) = T(x, v) 


Where 


x=0 

N-l 


T O, v) = ^ /(r, y) e ] 




y= o 
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Item 

Correct choice 

A 

IV 

B 

1 

C 

II 

D 

III 

E 

1 



S. The two ways are: 

• Forward mapping: for each pixel in the input image, find its location in the output image 
and assign its value 

o Multiple output values are assigned to the same output pixel 
o Some output locations may not assigned values. Hence, interpolation is necessary 
to fill these locations.' 

• Inverse mapping: for each pixel in the output image, find the location in the input image 
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on theinten^ir transform ' use the interpolation to calculate intensity value based 

on the mtensit.es in the input image location. 


I i . . ! «t>'- iwv-auuii. 

P efer doing inverse mapping is possible because the i 
input image pixel values. 


interpolation is directly calculated from the 


ti re S |S rat| ° n is the process of aligning different images taken for the same scene at different 
Th^h 53 * 6 lte ima ^ e ^ 0r mon ' tor ' n g earth movement) with different viewing angles/ distances. 

e o jective is to align the images to some reference images so that measurements and 
comparison could be made across all the images. 

A simple mode for registration is to build the mapping from a reference image f(x,y) to image 
g(v,w) using a set of tie points as follows: 

X — c ± v + c 2 w + C 3 VW + c 4 
y = c 5 v + c 6 w + c 7 vw -f c 8 

By Applying these two equations at four tie points to obtain 8 equations in eight unknowns, solve 
them for the unknowns, we get the model. The model is then used to transform all input image 
pixels. After transforming the coordinates, interpolation is used to assign intensity values 
The following figure explains image registration process. 
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